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/.  INTRODUCTION 


The  analysis  of  laser  speckle  interferograms  is  often  a  tedious  and  time  consuming  task. 
Many  attempts  have  been  made  toward  automating  the  analysis  process  with  excellent 
success.  The  only  major  drawbacks  have  been  in  the  time  required  and  cost  effectiveness  of  the 
process.  With  the  idea  of  reducing  cost,  decreasing  processing  time  and  increasing  reliability  a 
new  technique  was  sought  for  interpreting  speckle  interferograms.  The  new  method  utilizes  an 
off-axis  photosensor  element  coupled  to  a  mini-computer  to  achieve  these  goals. 

Laser  speckle  interferograms  are  most  commonly  used  to  make  deformation  measurements 
of  deformable  bodies.'  Figure  I  illustrates  the  basic  method  for  making  a  laser  speckle 
interferogram.  When  a  diffuse  surface  of  a  structure  is  illuminated  with  coherent  radiation,  a 
grainy  speckle  effect  is  imaged  by  the  eye  or  film  plane  of  a  camera  due  to  the  interference  of 
light  from  the  structure.  This  speckle  effect  is  enhanced  when  the  structure  has  microscopic 
surface  irregularities.  If  the  optical  configuration  remains  fixed,  the  speckle  pattern  of  the  test 
object  may  be  recorded  on  the  film  plane  of  a  camera.  Further,  if  the  structure  is  deformed,  the 
speckle  points  shift  with  the  deformation  and  a  second  exposure  of  the  deformed  speckle 
pattern  can  be  made. 

Using  a  technique  of  double  exposure,  speckle  interferograms  of  a  structure  are  normally 
made  by  photographing  the  speckle  pattern  in  a  deformed  and  undeformed  configuration.  A 
beam  of  laser  light  is  then  passed  through  a  region  of  the  double  exposure  where  the  local 
deformation  is  desired.  As  the  beam  passes  through  the  film,  the  deformed  and  undeformed 
speckle  recorded  there  diffract  the  laser  light  and  cause  an  interference  effect  on  a  viewing 
screen.  A  diffraction  halo  modulated  by  light  and  dark  bars  of  light  is  produced  where  the 
distance  2d  between  bars  is  inversely  proportional  to  the  distance  between  the  undeformed 
and  deformed  speckle  on  tlx  film  plane.  A  normal  to  the  light  and  dark  bar  pattern  indicates 
the  axis  of  deformation  of  the  speckle.  The  theory  to  be  developed  assumes  that  the 
deformation  region  illuminated  by  the  laser  beam  in  reconstruction  is  uniform  and  that  the 
linear  optical  theory  is  applicable.^ 

Although  the  laser  speckle  interferograms  can  be  made  rather  quickly,  the  analysis  time  of  a 
single  illuminated  point  can  be  very  significant,  i.e.,  several  seconds.  When  a  full  field  analysis 
is  required,  the  time  limitation  can  be  severe  if  several  thousand  points  of  deformation  data  are 

1.  J.  A.  SchaefTcl.  Acoustical  Speckle  Imer/eromelry,  VS  Army  Missile  Research  and  Development  Command.  Redstone 

Arsenal.  Alabama.  Technical  Report  T-79-39.  22  March  1979. 

2.  J.A.  Schaeffel,  B.R.  Mullinix,  W.F.  Ranson,  and  W.F.  Sy/imon.Compuier  Aided  Optical  Nondestructive  Flaw  Detection 

System  for  Composite  Materials,  US  Army  Missile  Research  and  Development  Command,  Redstone  Arsenal,  Alabama. 

Technical  Report  T-78-S,  26  September  1977. 
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needed.  Therefore,  a  rapid  low-cost  device  is  needed  which  can  generate  contours  of  constant 
displacement  over  a  full-field  interferogram. 

2.  THEORETICAL  DEVELOPMENT 

2  illustrates  the  reconstructed  diffraction  halo  modulated  by  light  and  dark  bars  of 
light.  From  the  linear  theory,^  the  displacement  in  the  6  direction  of  a  point  on  the  body  is 
given  as: 

SXf 

“9  “  2d  (I) 


Figur*  2.  Diffraction  halo  geomatry. 
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where. 


S  =  film  scale  factor  (magnification  ratio). 

A  =  wavelength  of  laser  illumination  source. 

f  =  distance  from  interferogram  to  analyzer  screen. 

d  =  distance  from  central  bright  spot  to  first  minima. 

displacement  of  the  point  illuminated  by  the  laser  on  the  object  in  the  0 
direction. 


The  vertical,  U,  and  horizontal,  Uh  components  of  displacement  may  be  obtained  from  U0  as:’ 
Ufl  “  Ug  CO80  ■  CO80  (2) 

U ,  -  U„  8in0  -  4^  8in0  (3) 

V  V  2d 

and  from  the  geometry. 


■  CO80 

-  8ln0 

therefore. 


SAf 

“  ^ 


SAf 


(4) 

(5) 


(6) 

(7) 


The  optical  configuration  of  Figure  3  shows  a  photodetector  placed  off-axis  a  distance  X 
from  the  laser  central  bright  spot  and  in  the  field  of  the  diffraction  halo.  The  photodetector 
senses  the  light  intensity.  Whenever  a  minima  in  intensity  is  detected,  the  fringe  order  at  the 
illumination  point  is  an  integral  number  n  =  1,  2,  3, ...  in  value. 


i.  J.A.  Schaeffel.  Ijib  Notebook  7232,  US  Army  Missile  Research  and  Development  Command,  Redstone  Arsenal. 
Alabama.  10  November  1978. 


By  moving  the  interferogram  in  a  raster  fashion  and  recording  the  intensity  at  each  point, 
contours  of  constant  n  are  generated.  The  first  minima  corresponds  to  n  =  I,  the  second 
minima  n  =  2,  etc.  These  contours  represent  loci  of  constant  displacement  on  the  body  from 
which  the  interferogram  was  made.  1  he  direction  of  displacement  is  parallel  to  the  axis  formed 
from  the  central  bright  spot  and  the  central  axis  of  the  photodetector.  Once  n  is  known  at  a 
particular  point,  the  displacement  along  the  direction  parallel  to  the  X  axis  is  given  as: 


where. 


V 


SU 

2X 


(9) 


an  optical  geometrical  constant. 

In  operation  the  interferogram  is  pointwise  scanned  back  and  forth  in  front  of  the 
photodetector.  For  each  point  of  interest  on  the  interferogram  the  intensity  is  recorded. 
Contours  of  minima  and  maxima  intensity  are  determined  from  this  record.  These  contours 
correspond  to  lines  of  constant  displacement  along  the  X  optical  axis  as  determined  by 
Equation  8. 

To  illustrate  the  method  of  analysis,  consider  a  uniaxial  tension  specimen  as  shown  in 
Figure  4.  The  specimen  is  loaded  in  the  x  direction  with  a  load  P  which  causes  a  deformation 
parallel  to  the  x  axis  neglecting  poisson  effects.  Suppose  an  interferogram  is  made  of  the 
structure  before  and  after  deformation.  The  interferogram  is  analyzed  with  a  photodetector 
offset  a  distance  X  in  the  x  direction  of  the  interferogram  and  the  intensity  is  shown  plotted  as  a 
function  of  x  and  Y  for  24  locations  on  the  tension  specimen.  At  each  location  on  the 
interferogram  intensity  plot  which  is  a  minimum,  the  displacement  is  an  integral  number 
n  =  0,  I,  2,  ...  of  p.  At  locations  where  the  intensity  is  a  maximum  the  displacement  is 
n  =  1/2,  3/2,  5/2, ...  of  v. 

In  a  practical  situation  the  intensity  reading  of  the  photodetector  is  digitized  and  stored  in 
digital  form  by  a  computer  for  maxima  and  minima  contour  determinations.  Due  to  the 
inherent  noise  found  in  the  laser  speckle  diffraction  halo,  some  digital  filtering  of  the  optical 
signal  is  desired.  This  may  be  accomplished  by  a  discrete  fourier  transform  (DFT)  of  the 
intensity  field.  Given  N  samples  of  the  intensity  field,  which  are  uniformly  separated  and 
collinear  along  an  axis  parallel  to  x  or  y.  shown  in  Figure  4,  when  the  frequency  domain  of  the 
transformed  field  is  given  as:^ 


4.  J.A.  Bells,  Signal  Processing  Modulation  and  Noise,  American  Elsevier  Publishing  Company,  Inc.,  New  York,  1971. 


Cn 


,N-1  -J  n2-,TK 

I(k)e  N 


(10) 


where  n  =  0,  1,  2, ...  N-l 


In  speckle  filtering  a  lowpass  harmonic  filter  is  desired  with  a  harmonic  cut-off  of  P  where  /3 
<  N-l.  The  larger  values  of  result  in  less  filtering  of  the  speckle  noise. 


The  final  filtered  intensity  field  I'  is  then  given  as: 

I'  (k)  -  Z  -  Cn  e 
n«a 


(II) 


Practice  dictates  that  jS  should  be  slightly  larger  than  the  maximum  number  of  contour 
minima  or  maxima  for  the  best  performance  of  the  digital  filter.  The  lower  harmonic  retained 
in  the  filter  is  determined  by  the  value  of  a  which  is  normally  0. 


3.  EXPERIMENTAL  CONFIGURATION 


The  mechanical  and  electrical  block  diagram  for  implementing  the  theory  of  Section  2  is 
shown  in  Figure  5,  and  6.  Figure  7  shows  a  photograph  of  the  complete  system.  The  basic 
components  and  their  functions  are: 

•  PDF  1 1  /40  Minicomputer  -  This  computer,  manufactured  by  the  Digital  Equipment 
Corporation,  (DEC),  serves  as  the  central  controller  for  analyzing  an  interferogram.  It 
contains  32,000  -  16  bit  words  of  storage  and  uses  a  DEC  RTl  1  software  control  system. 

•  RXOl  Floppy  Disk  -  The  disk  drive  is  used  to  store  programs,  data  and  the  RTIl 
monitor.  It  is  a  dual  drive  unit  with  250,000-16  bit  words  of  storage.  It  is  manufactured  by 
DEC. 


•  DECWRITER  -  This  is  the  typewriter  input-output  port  for  the  computer.  It  is  a  30 
character  per  second  unit  used  to  input  scan  parameters  for  analyzing  an  interferogram  and 
printing  the  results  of  the  scanning  process.  It  is  manufactured  by  DEC. 

•  ARM  -  This  unit  is  the  analog  input  port  manufactured  by  DEC.  It  is  used  to  digitize  the 
photodetector  intensity  signal.  The  unit  features  a  16-channel  A/D  converter  with  an  input 
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Figure  6.  Experimental  electrical  configuration. 
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voltage  range  of  0  to  5  V,  a  resolution  of  10  bits,  and  an  accuracy  of  ±0. 1%  of  full  scale.  It 
also  has  a  22  to  24  ^sec  conversion  time,  an  input  impedance  of  10  Mfl  and  a  maximum 
settling  time  of  iftscc. 

•  Oscilloscope  A  -  This  is  a  Tektronix  Company  Model  549  oscilloscope  with  a  type  W 
differential  comparator  used  as  a  high  gain  amplifier  to  amplify  the  photodetector  signal. 

•  Oscilloscope  B  -  This  is  a  Hewlett-Packard  Model  HPI700B  oscilloscope  used  to 
monitor  the  input  voltage  to  the  ARI I  A/D  converter. 

•  Photodetector  -  The  photodetector  selected  is  a  Spectra  Physics  Model  40IC 
photodetector  with  a  peak  spectral  sensitivity  at  6328  angstroms. 

•  LowPass  Filter/  Voltage  Clipper  -  This  circuit  consists  of  a  lowpass  frequency  passive 
filter  and  voltage  clipper.  The  RC  time  constant  is  0.1 445  OF  with  a  cut-off  clip  voltage  of 
about  3.8  V.  The  circuit  is  used  to  reduce  signal  noise  and  prevent  a  voltage  overload  condition 
to  the  ARI  1. 

•  DR  I  IC-The  DR  1 1C  is  a  digital  I/O  interface  used  to  position  xY-table  for  scanning 
interferograms.  It  is  manufactured  by  DEC  and  features  16  ports  of  input  and  16  ports  of 
output  TTL  compatible  voltages. 

•  XY  -  Table  Controllers  -  Control  of  the  XY-Table  via  the  DRl  1C  interface  is 
maintained  by  these  units.  This  system  contains  two  slo-syn  translator  modules  (Type 
STM  101)  manufactured  by  the  Superior  Electric  Company.  Also  included  in  the  system  are 
two  type  M  PS- 1 000  power  supplies  also  made  by  Superior. 

•  XY-table  -  The  XY-table  used  to  position  the  interferogram  is  made  by  Design 
Components,  Inc.,  and  is  a  Model  DC-66  having  a  scanning  range  of  152.4  mm  in  the  x  and  Y 
directions.  It  can  take  single  steps  between  scan  points  of  0.0254  mm  with  a  repositioning 
accuracy  of  0.00254  mm. 

•  He-Ne  laser  -  The  laser  selected  was  a  Spectra-Physics  Model  124  He-Ne  laser  operating 
at  6328  A. 

Operation  of  the  optical  system  is  controlled  through  the  PDP  1 1/40  computer.  A 
description  and  listing  of  the  computer  code  is  contained  in  the  appendix.  After  the  computer 
program  is  started,  the  parameters  for  performing  a  raster  scan  of  the  interferogram  are  input 
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to  the  computer  through  the  DEC'WRI TFR.  A  total  of  KM)  points  in  the  x  direction  and  50 
points  in  the  y  direction  may  be  scanned  for  a  total  of  5000  points.  After  the  scanning 
parameters  are  input  to  the  computer,  a  request  is  made  for  the  initial  position  for  scanning. 
The  operator  then  positions  the  interferogram  interactively  with  the  computer.  Once  this 
operation  is  complete  the  interferogram  is  scanned  by  moving  it  to  each  scan  point 
automatically  and  the  photodetector  intensity  signal  is  digitized  by  the  computer.  This  entire 
operation  is  automatic.  If,  during  the  scanning  process,  the  operator  needs  to  stop  scanning 
and  return  to  the  initial  scan  position,  the  setting  of  any  bit  on  the  control  console  may 
perform  this  function.  After  the  scan  is  made,  a  DFT  filter  operation  is  performed  if  requested, 
and  a  contour  plot  of  the  results  is  made.  While  the  contour  plot  is  being  made,  the  results  can 
be  modified  by  selecting  the  proper  bit  configuration  on  the  computer  console.  A  listing  of 
these  commands  is  contained  in  the  appendix. 

The  method  of  contour  plotting  selected  involves  the  use  of  a  variable  reference  signal 
intensity  plane.  Intensity  information  from  the  photodetector  is  stored  as  a  digital  signal 
floating  point  level  from  0.0  to  1.0  in  the  computer  memory  and  may  be  represented  in 
graphical  form  as  in  Figure  4.  To  make  a  contour  plot  of  deformation,  the  operator  specifies  a 
signal  ratio  from  0.0  to  1 .0.  All  intensities  stored  in  memory  greater  than  the  signal  ratio  are 
plotted  as  light  regions  (maxima)  and  all  intensities  less  than  the  signal  ratio  are  plotted  as 
dark  regions. 

Smoothing  of  the  contours  is  possible  using  the  DFT  filters  provided.  The  filter  may 
precede  or  follow  after  a  contour  plot  is  made. 

4.  EXPERIMENTAL  RESULTS 

In  order  to  test  the  operation  of  the  system  and  new  scanning  technique,  a  composite  E-glass 
tube  was  fabricated  using  an  Epon  828  resin  system.  Figure  8  shows  a  schematic  of  the  tube. 
Pertinent  tube  data  includes  the  following; 

•  Cylinder  No.  8-E-45-A  •  74. 1 1  mm  inside  diameter 

•  Test  Pressures  0-150  psig  •  1.27  mm  wall  thickness 

•  Cut  Fiber  Flaw  •  76.65  mm  outside  diameter 

•  E-glass,  Epon  828  resin  system  •  75.38  mm  mean  diameter 
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The  tube  had  an  overall  length  of  2S3.49  mm  with  a  cut-flaw  in  the  fiber  126.75  mm  from  the 
bottom  of  the  tube.  1  he  tube  was  painted  in  the  region  marked  E  for  a  total  of  202.95  mm. 
Specimen  pressure  testing  was  at  0  and  1 50  psig.  Boundary  conditions  on  the  endcaps  included 
a  fixed-end  condition  at  the  base  and  a  free-floating-end  condition  at  the  top  of  the  tube.  After 
pressurizing  the  tube  to  1 50  psig,  a  laser  photograph  was  made  in  the  deformed  condition  and 
then  the  pressure  was  released.  A  second  laser  photograph  was  then  made  for  the  double 
exposure. 

Analysis  was  done  with  a  total  ot  100  points  in  the  x  direction  and  50  points  in  the  Y 
direction.  The  raster  scan  included  a  total  of  5000  points.  For  the  scanning  process  X  =  152.4 
mm,  f=255.27  cm  and  v=0.0l209  mm/  fringe  order.  The  interferogram  scanning  increments  in 
X  and  y  covered  the  entire  painted  region  uniformly.  Inteferograms  of  the  specimen  were 
analyzed  for  the  front  and  back  surfaces  with  and  without  a  flaw,  respectively. 

Figures  9  and  10  are  contour  displacement  plots  of  the  flawed  cylinder  surface  made  with 
the  scanning  system.  Maxima  are  indicated  with  a  (-)  and  minima  with  an  (*).  The  horizontal 
row  of  numbers  indicates  the  x  axis  scan  position  and  the  vertical  column  indicates  the  y  axis 
scan  position.  The  flawed  region  of  the  tube  is  indicated  by  a  broad  minima  at  location  column 
60,  row  25.  Figures  II  and  12  indicate  the  results  after  using  a  DFT  Filter.  The  filter 
parameters  are  o  =  0  and  =  14  in  both  x  and  y  directions.  As  seen  in  Figures  II  and  12, 
there  is  a  significant  perturbation  in  the  displacement  contour  map  at  the  same  location  as  in 
Figures  9  and  10.  The  results  of  the  filter  are  to  reduce  the  overall  signal  noise  level.  The 
difference  in  displacement  between  a  maxima  and  a  minima  contour  is  0.00605  mm. 
Therefore,  high  strain  gradients  are  present  in  the  tube  at  those  locations  where  displacement 
contours  have  sharp  changes  in  the  x  or  y  directions. 

A  series  of  scans  followed  by  manual  fringe  digitizations  were  made  and  the  results 
compared.  In  general,  the  percentage  of  difference  between  the  automated  and  manual 
analysis  differed  by  less  than  4%  for  fringe  orders  of  approximately  n  »  1.  As  n  approaches 
values  of  4  to  5  the  measured  displacements  differ  by  less  than  2%.  The  system  is  very  fast  and 
has  a  potential  digitization  rate  of  6250  data  points/ min.  Under  certain  optimum  conditions 
10*  data  points  may  be  digitized  in  about  2  min  if  the  distances  between  scan  points  is 
minimized. 

An  interesting  test  was  made  by  plotting  on  an  oscilloscope  the  intensity  map  as  digitized  by 
the  computer  and  without  filtering.  Figure  13  shows  the  result  which  may  be  compared  with 
Figures  9  and  70.  In  this  test,  a  composite  tube  with  a  similar  flaw  as  before  was  analyzed  and 
discontinuities  in  the  fringe  pattern  can  be  observed. 
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Figure  9.  Rawed  cylinder  displacement  contour  map  without  discrete  Fourier 
Transform  Filter  (front  view).  Signal  ratio  =  0.5. 
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Flawed  cylinder  displacement  contour  map  with 
Transform  Filter  (front  view).  Signal  ratio  =  0.6. 
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Figure  11.  Flawed  cylinder  displacement  contour  map  with  discrete  Fourier 
Transform  Filter  (front  view).  Signal  ratio  =  0.5. 
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Figure  12.  Flawed  cylinder  displacement  contour  map  with  discrete  Fourier 
Transform  Filter  (front  view).  Signal  ratio  =  0.6. 


Figure  13.  Displacement  contour  mapping  of  a  flawed  composite  cylinder. 
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5.  CONCLUSIONS 


The  theory  and  an  experimental  configuration  were  developed  for  automated  displacement 
contour  analysis  of  laser  speckle  interferograms.  The  number  of  displacement  contours 
available  is  proportional  to  the  number  of  diffraction  fringe  orders  n  available  throughout  the 
offset  distance  X.  The  system  was  checked  to  determine  how  well  fringe  contours  could  be 
measured  and  found  to  be  accurate  to  at  least  4%  of  the  measured  distance  by  hand.  System 
scanning  speed  is  extremely  important,  and  it  was  found  that  5000  data  points  could  be 
analyzed  in  about  4  min  when  the  system  was  operating  at  about  20%  of  its  maximum  rated 
speed.  This  is  very  encouraging  for  systems  in  which  as  many  as  10*  data  points  would  be 
analyzed.  Advantages  of  the  present  system  over  those  of  the  past  include  virtual  operator-free 
scanning,  high  speed,  low  cost  and  full  field  displacement  field  mapping. 


APPENDIX 


COMPUTER  INPUT  VARIABLES 


I.  CONSOLE  CONTROL  SWITCHES  (POP  11/40): 


Bit  0  =  A  new  signal  ratio  will  be  requested  for  making  a  new  contour  plot  immediately  if 
this  switch  is  turned  on  while  a  contour  plot  is  being  made. 

Bit  1  =  A  new  signal  ratio  will  be  requested  for  making  a  new  contour  plot  at  the  end  of 
the  current  contour  plot  if  this  switch  is  turned  on. 

Bit  2  =  A  new  Discrete  Fourier  Transform  filter  parameter  for  filtering  the  current 
intensity  signal  will  be  requested  at  the  end  of  the  current  contour  plot  if  this 
switch  is  turned  on. 

Bit  3  =  If  turned  on  while  contour  plotting,  a  new  set  of  scan  parameters  will  be  requested 
at  the  end  of  the  current  contour  plot. 

Bit  4  =  If  turned  on  while  contour  plotting,  a  set  of  new  ASCII  code  plot  symbols  will  be 
requested  at  the  end  of  the  current  contour  plot. 

Bit  S  =  If  turned  on,  fringe  orders  may  be  ordered  to  displacement  contours  at  the  end  of 
the  current  contour  plot. 

Bit  6  =  End  the  program  if  turned  on  while  contour  plotting.  The  current  contour  plot  is 
completed. 

2.  CODE  VARIABLES. 

IXP  =  total  number  of  points  in  the  x  direction  to  be  scanned. 

IXS  =  number  of  0.02S4  mm  increments  between  scan  points  in  the  x  direction. 

lYP  =  total  number  of  points  in  the  y  direction  to  be  scanned. 

lYS  =  number  of  0.0254  mm  increments  between  scan  points  in  the  y  direction. 

IXFI  =0+1  for  the  x  direction. 

IXF2  =^+l  for  the  x  direction. 
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lYFI  =  a+1  for  the  y  direction. 


1YF2  =  ^+1  for  the  y  direction. 

SR  =  signal  ratio. 

1 1  =  number  of  0.0254  mm  increments  to  advance  the  x  stage  before  scanning. 

12  =  number  of  0.0254  mm  increments  to  advance  the  y  stage  before  scanning. 


C - PROGRAM  SOFT 

C - WRITTEN  BY5  JOHN  A.  SCHAEFFELr  JR. 

C - BIT  0«NEU  SR  AT  END  OF  PRINT  LINE 

C - BIT  1=NEU  SR  AT  END  OF  CONTOUR  PLOT 

C - BIT  2=NEU  DFT  AT  END  OF  CONTOUR  PLOT 

C - BIT  3«NEU  SCAN  PARAMETERS 

C - BIT  4=NEW  ASCII  PLOT  CODES 

C - BIT  5*EXAMINE  FRINGES 

C - BIT  6=END  PROGRAM 

C - ALL  BITS  MUST  BE  OFF  FOR  SCAN 

C - 0-5  VOLTS  ARll  RANGE 

C - INPUT  STREAM  ON  AR11>3 

DIMENSION  A(100>S0>>B(2f 100>rIT(101) 

KA1=72 
KA2=76 
URITE(5f 1) 

1  FORMAT!'  SDFT  ANALYZER  CODE') 

41  URITE(Sf2> 

2  FORMAT!'  NO.  X  POINTS.  X-SCAN  INTERVAL-213' > 

READ!S.3)  IXP.IXS 

3  F0RMAT!213> 

URITE!S>4) 

4  FORMAT!'  NO.  Y  POINTS.  Y-SCAN  INTERVAL-213') 

READ!S.3)  lYP.IYS 

URITE!S.ll) 

11  FORMAT!'  DFT  X.Y  MIN-MAX  FILTER  LEVELS-4I3') 

REA0!5.12)  IXF1.IXF2.IYF1.IYF2 

12  F0RMAT!4I3) 

URITE!S.13) 

13  FORMAT!'  SIGNAL  RATIO-FIO.O' ) 

READ!5.14)  SR 

14  F0RMAT!F10.0) 

10  URITE!S.36) 

36  FORMAT!'  X.Y  STAGE  POSITION-214') 

READ!5.37)  11.12 

37  FORMAT! 214) 

IF!11.EO.O.AND.12.EQ.O)  GOTO  38 
CALL  XA0V!I1.S) 

CALL  YADV!I2.S) 

GOTO  10 

38  CONTINUE 
IM-0 
IYl-1 
IY2-IYP 
JY3-1 

IXS-IAB8!IXS) 

IYS-IAB8!IYS) 

DO  6  I-l.IXP.l 

DO  5  J-IY1.IY2.IY3 

CALL  IPOKE!* 170400. *21001) 

ITP-IPEEK!* 170402) 

A!  I  .  J)=.FLOAT!  ITP)/1024. 

1TEST«IPEEK! *177570) 

IF!ITEST.NE.O)  GOTO  41 
IF!J.NE.IY2)  CALL  YADV!IY8.S) 

5  CONTINUE 
IY4«IY1 
IV1-IY2 
IY2-IY4 
IY3— IY3 
IY8— IY8 

CALL  XAOV!IXS.S) 

6  CONTINUE 
1X8— 1X8 
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DO  7 

CALL  XADVdXStS) 

7  CONTINUE 
irdYS.GT.O)  GOTO  9 
DO  e  J=ltIYP-l.l 
CALL  YADVdYSfS) 

8  CONTINUE 

9  CONTINUE 
IFdXFl.EQ.O)  GOTO  39 
IN-0 

C-6.2B318S 

D=FLOATdXP) 

DO  17  J=ltIYPrl 
DO  IS  I-ldXP>l 
Bdf  t)<0. 

15  B(2>I)-0. 

DO  14  K-IXF1-1»IXF2-1»1 
DO  16  L-0»IXP-lfl 

B<ltK-fl)-B(lrK-fl>'f(A(L+l>J)«COS(FLOAT(K»L>«C/D)  >/D 
B(2>K'fl>-B(2»K+l>-(A<LFlf  J>«SIN(FLOAT<K«L>«C/D)  >/D 

16  CONTINUE 

DO  17  H-OdXP-ld 
A(N'flf  J>-0. 

DO  17  N-IXFl-lf IXF2-1f1 

A<M'df  J>-A(N-df  J>+B(lfN-t-l>«COS(FLOAT<N«N)«C/D> 
A(N'flFj)>A(N'fl>J)-B<2>N+l>«SIN<FL0AT(N«N)«C/D) 

17  CONTINUE 
D»FLOATdYP) 

DO  20  J-lrlXPrl 
DO  18  I-l.IYPfl 
BdrI)-0. 

18  B(2fI>-0. 

DO  19  K-IYFl-l»iyF2-lf 1 
DO  19  L-0»IYP-1.1 

B(lrK-fl>-B(l>KFl>  +  <A(JrL'fl)«COS<FLOAT(K»L)«C/D>)/D 
B(2fK+l)-B(2FK-fl>-(A<JFLFl>«8IN(FL0AT(K»L>«C/D>)/D 

19  CONTINUE 

DO  20  M-0.IYP-1»1 
A<J>N'fl)-0. 

DO  20  N-IYF1-I»IYF2-1»1 

A(J>H-»'l>-A(jFM-fl>'FB(l>NFl)«COS(FLOAT(N«N>»C/D> 

A<  JFH'«-l)-A(jFN'l-l)-B<2rNfl>»8IN<FL0AT<N«N>«C/D) 

20  CONTINUE 

39  CONTINUE 
IFdN.EO.O)  GOTO  40 
WRITE(Sfl3> 

READ(Sfl4)  SR 

40  CONTINUE 
ANIN-A<1>1) 

ANAX-Adrl) 

DO  21  I-1>IXPf1 
DO  21  J-lfIYP>l 

IF(A(IfJ).OT.AHAX)  AHAX-AdfJ) 

21  IF(A(IfJ)  .LT.AHIN)  ANIN-Ad>J) 

IF ( ANAX . EQ . AHIN )  ANAX-ANAXt 1 . 

URITE(5>22) 

22  FORHAT('  DI8PLACENENT  FIELD  NAPt') 

DO  23  I-lflOfl 

DO  23  J-1>10>1 
22  IT<(I-l)«104J>-J'f47 

IT<101)-48 

URITE(5f24>  (ITlJltJ-lflXP-fltl) 

24  FORMAT (IN  fIOIAI) 

DO  25  I-l>IXPfl 

25  IT(I>-4A 
KCR-49 
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URITE(5.?6>  KCR.<IT( J> • J=lf IXPf 1> 

26  FORMATdH  tlAlrlOOAl) 

VO  33  I=IYP-1.2.-1 

KCR=KCR+1 

no  31  M=2.IXP-1.1 

31  IT(M)=46 

no  32  J=2.IXP-1.1 
Al=>(A(J>I)-Af1IN)/(AHAX-AMIN) 

IFCAl.GT.SR)  IT(J)»KA1 

32  IFXAl.LE.SR)  II(J)>KA2 

WRITE! 5. 26)  KCR f  ( IT ( J)  .  J«^l ,  IXP .  1 ) 

IF<KCR.EQ.57)  KCR«47 

IM=0 

ITEST=IPE£K< ’177570) 

IF(ITEST.NE.O)  IM=1 
IFdTEST.EQ.l)  GOTO  39 

33  CONTINUE 

DO  34  I^lrlXPfl 

34  IT(I)c46 
KCR-KCR+1 

WRITE (5> 26)  KCRf  (IT(J)  IXPtl) 

IF<ITEST.EQ.2)  GOTO  39 

IM-0 

IF(ITEST.EQ.4)  WRITE(Sdl) 

IF(ITEST.EQ.4)  READ(Sfl2)  IXFl f IXF2f I YFl r I YF2 
IF(ITEST.E0.4)  GOTO  9 
IFdTEST.EQ.S)  GOTO  41 
IF<ITEST.EQ.16)  WRITE(S>42> 

42  FORMAT!'  HIGHr  LOW  ASCII  CODeS-212') 

IF!ITEST.E0.16)  READ!5f43)  KA1>KA2 

43  FORMAT! 212) 

IF!ITEST.Ea.l6)  GOTO  40 
IF!ITEST.CQ.64)  GOTO  35 
IF!ITEST.NE.32)  GOTO  10 
1BX»1 

IBY=IYP 

IX8«IABS!IXS> 

IYS»IABS!IYS> 

44  WRITE!5t4S) 

45  FORMAT!'  X»Y  INPUT  POSITION-213') 

READ!5t3)  IZXdZY 
I1'°!IZX-IBX)«IXS 

12— !IZY-IBY>*iyS 
IF!Il.NE.O)  CALL  XADV!I1>5> 

IF!I2.NE.O)  CALL  YADV!12>5) 

IBX-IZX 

IBY»IZY 

IFdZX.NE.l.OR.IZY.NE.iyP)  GOTO  44 
GOTO  10 

35  CONTINUE 
STOP 
END 

SUBROUTINE  XADU!ISfIR> 

X»0. 

IFdS.OT.O)  GOTO  3 
IP-1ABS(IS> 

DO  2  I-l>IPfl 

CALL  IP0KE!*167772»’100000> 

DO  7  K-lflRrl 
7  Y-8IN!X> 

CALL  IP0K£!*167772> ’000000) 

DO  1  J>lfIR>l 

1  Y-8IN!X) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 
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DO  5  11=1. IS. 1 
CALL  IPOKEC 167772, ’OAOOOO} 
DO  8  KK=l.IR.l 
8  y=SIN(X) 

CALL  IPOKE( *167772. '000000) 
DO  4  JJ=l.IR.l 

4  r=SIN<X) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 

SUBROUTINE  TADVaS.IR) 

X=0. 

IFCIS.GT.O)  GOTO  3 
IP=IABS(IS) 

DO  2  1=1. IP. 1 

CALL  IPOKEf *167772. *020000) 
DO  7  K=l.IR.l 

7  Y»SIN(X) 

CALL  IPOKE< *167772. *000000) 
DO  1  J'lfIR.l 

1  Y=SIN<X) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  II=l>IS.l 

CALL  IPOKE< *167772. *010000) 

DO  8  KK=1>1R.1 

8  Y»SIN(X) 

CALL  IPOKE<* 167772. *000000) 
DO  4  JJ=l.IR.l 

4  Y=SIN<X) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 

« 
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LIST  OF  SYMBOLS 


Symbol 

Cn 

2d 

dv 

dH 

f 

I(K) 

r(K) 

J 

N 

n 

S 

U« 

Uv 

UH 


Definition 

Discrete  Fourier  transform  frequency  spectral  components 

Spacing  between  n  =  1/2  fringe  orders 

Vertical  vectoral  component  of  d 

Horizontal  vectoral  component  of  d 

Distance  from  interferogram  to  analyzer  screen 

Digitized  intensity  data 

Filtered  intensity  data 

Imaginary  number 

Total  number  of  data  samples 

Fringe  order  of  contour 

Film  scale  factor 

Absolute  displacement  along  6  axis 
Absolute  displacement  along  vertical  axis 
Absolute  displacement  along  horizontal  axis 
Lower  harmonic  level  of  DFT  filter 
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Upper  harmonic  level  for  DFT  filter 


P 

X 

e 


Wavelength  of  laser  light  source 
Angle 

Optical  geometrical  constant 
Photodetector  offset  distance 
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